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Complex spectra of random matrices are studied in the regime of weak non- 
Hermiticity. The matrices we consider are of the form LLi + iH2, where Hi and LL2 
are Hermitian and statistically independent. In the first part of the paper we con- 
sider the case of matrices LL\^2 having i.i.d. entries. For such matrices the regime of 
weak non-Hermiticity is defined in the limit of large matrix dimension A'' by the con- 
dition {TyHI) (X 7V(Tr//|>. We show that in the regime of weak non-Hermiticity the 
distribution of complex eigenvalues of Hi+iH2 is dictated by the global symmetries 
of L[i,2, but otherwise is universal, i.e. independent of the particular distributions 
of their entries. Our heuristic proof is based on the supersymmetric technique and 
extends also to "invariant" ensembles of i/i,2. In the second part of the paper we 
study Gaussian complex matrices in the regime of weak non-Hermiticity. Using the 
mathematically rigorous method of orthogonal polynomials we find the eigenvalue 
correlation functions. This allows us to obtain explicitly various eigenvalue statis- 
tics. These statistics describe a crossover from Hermitian matrices characterized by 
the Wigner-Dyson statistics of real eigenvalues to strongly non-Hermitian matrices 
whose complex eigenvalues were studied by Ginibre. Two-point statistical measures 
such as spectral form-factor, number variance and small distance behavior of the 
nearest neighbor distance distribution p{s) are studied thoroughly. In particular, 
we found that the latter function may exhibit unusual behavior p{s) oc s^^^ for some 
parameter values. 



I. INTRODUCTION 

Eigenvalues of large random matrices have been attracting much interest in theoretical physics 
since the 1950's [0-0. Until recently only the real eigenvalues were seen as physically relevant, 
hence most of the studies ignored matrices with complex eigenvalues. Powerful techniques to deal 
with real eigenvalues were developed and their statistical properties are well understood nowadays 
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0. Microscopic justifications of the use of random matrices for describing the universal spectral 
properties of quantum chaotic systems have been provided by several groups recently, based both on 
traditional semiclassical periodic orbit expansions |^,^ and on advanced field-theoretical methods 
pO| , pl| . These facts make the theory of random Hermitian matrices a powerful and versatile tool 
of research in different branches of modern theoretical physics, see e.g. 



Recent studies of dissipative quantum maps ^2 13 1, asymmetric neural networks |14 1^, and 
open quantum systems [^-po[ stimulated interest to complex eigenvalues of random matrices. 
Most obvious motivation comes from studies of resonances in open quantum systems, i.e. systems 
whose fragments can escape to or come from infinity. The resonances are determined as poles 
of the scattering matrix (S- matrix), as a function of energy of incoming waves, in the complex 
energy plane. The real part of the pole is the resonance energy and the imaginary part is the 
resonance half-width. Finite width implies finite life-time of the corresponding states. In the 
chaotic regime the resonances are dense and placed irregularly in the complex plane. Recently, the 
progress in numerical techniques and computational facilities made available resonance patterns of 
high accuracy for realistic open quantum chaotic systems like atoms and molecules 

Due to irregularity in the resonance widths and positions the S-matrix shows irregular fluc- 
tuations with energy and the main goal of the theory of the chaotic scattering is to provide an 
adequate statistical description of such a behavior. The so-called "Heidelberg approach" to this 



problem suggested in 1 22 1 makes use of random matrices. The starting point is a representation of 
the S-matrix in terms of an effective non-Hermitian Hamiltonian Ti-eff — H — iT. The Hermitian 
N X N matrix H describes the closed counterpart of the open system and the skew-Hermitian 
iT — iWW^ arises due to coupling to open scattering channels a = 1, . . . , AT, the matrix elements 
Wja being the amplitudes of direct transitions from "internal" states i = 1, . . . , to one of open 
channels. The poles of the S-matrix coincide with the eigenvalues of Ti-eff - In the chaotic regime 
one replaces H with an ensemble of random matrices of an appropriate symmetry. This step is 
usually "justified" by the common belief according to which the universal features of the chaotic 
quantum systems survive such a replacement As a result, various features of chaotic quan- 

tum scattering can be efficiently studied by performing the ensemble averaging. The approach has 
proved to be very fruitful (for an account of recent developments see [po[). In particular, it allowed 
to obtain explicitly the distribution of the resonances in the complex plane for chaotic quantum 
systems with broken time-reversal invariance [^,^ and, in its turn, this distribution was used to 
clarify some aspects of the relaxation processes in quantum chaotic systems p3| . 

A very recent outburst of interest to the non-Hermitian problems p^|3^ deserves to be men- 
tioned separately. During the last several years complex spectra of random matrices and operators 
emerged in a diversity of problems. Hatano and Nelson described depinning of flux lines from 
columnar defects in superconductors in terms of a localization-delocalization transition in non- 
Hermitian quantum mechanics . Their work motivated a series of studies of the corresponding 
non-Hermitian Schrodinger operator p7|-|31| and, surprisingly, random matrices appeared to be 



relevant in this context [g^JcH • Complex eigenvalues were also discussed in the context of lattice 
QCD. The lattice Dirac operator entering the QCD partition function is non-Hcrmitian at nonzero 
chemical potential and proves to be difficult to deal with both numerically and analytically. Recent 
studies of chiral symmetry breaking used a non-Hermitian random matrix substitute for the Dirac 
operator . There exist also interesting links between complex eigenvalues of random matrices 
and systems of interacting particles in one and two spatial dimensions ]3^ . And, finally, we have 
to mention that random matrices can be used for visualization of the pseudospectra of non-random 
convection-diffusion operators ||3^, and for description of two- level systems coupled to the noise 
reservoir js^. 

Traditional mathematical treatment of random matrices with no symmetry conditions imposed 
goes back to the pioneering work by Ginibre ||3^ who determined all the eigenvalue correlation 
functions in an ensemble of complex matrices with Gaussian entries. The progress in the field 
was rather slow but steady P,^7|-^, see also In addition to the traditional approach 

other aproaches have been developed and tested on new classes of non-Hermitian random matrices 



| [15| , |17 19 44|-[48|. However, our knowledge of the statistical properties of complex eigenvalues of 
random matrices is still far from being complete, in particular little is known about the universality 
classes of the obtained eigenvalue statistics. 

When speaking about universality one has to specify the energy scale, for the degree of univer- 
sality depends usually upon the chosen scale. There exist two characteristic scales in the random 
matrix spectra: the global one and the local one. The global scale is aimed at description of 
the distribution of the eigenvalues in bulk. The local one is aimed at decription of the statistical 
properties of small eigenvalue sets. For real spectra, the global scale is that on which a spectral 
interval of unit length contains on average a large, proportional to the matrix dimension, number 
of eigenvalues. If the spectrum is supported in a finite interval [a, b] the global scale is simply given 
by the length of this interval. In contrary, the local scale is that determined by the mean distance 
A between two neighbouring eigenvalues. Loosely speaking, the local scale is N times smaller than 
the global one sufficiently far from the spectrum edges, N being the matrix dimension. 

Universality in the real spectra is well established. The global scale universality is specific to 
random matrices with independent entries and does not extend to other classes of random matrices. 
The best known example of such universality is provided by the Wigner semicircle law fist] : 

= ^v/4J2-X2 = Ni^iX) = i, (1) 

which holds for random matrices whose entries satisfy a Lindeberg type condition | |50[ |. In this 
expression the parameter J just sets the global scale in a sense as defined above. It is determined 
by the expectation value — Tri?^). It is generally accepted to scale entries in such a way 
that J stays finite when N oo, the local spacing between eigenvalues in the neighbourhood of 
the point X being therefore A cx Similar universality is also known for complex spectra 

p2|,p|. 
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From the point of view of universality the semicircular eigenvalue density is not extremely robust. 
Most easily one violates it by considering an important class of so-called "invariant ensembles" 
characterized by a probability density of the form ViH) cx exp {— N Tr V (H)) , with V{H) being 
an even polynomial. The corresponding eigenvalue density turns out to be highly nonuniversal 
and determined by the particular form of the potential V{H) ||5l|,|5^. Only for V{H) = it is 
given by the semicircular law, Eq.(^. Moreover, one can easily have a non-semicircular eigenvalue 
density even for real symmetric matrices S; Sij — Sji with i.i.d. entries, if one keeps the mean 
number of non-zero entries p per column to be of the order of unity when performing the limit 
N — !■ oo. This is a characteristic feature of the so-called sparse random matrices [^3|-^. 

Much more profound universality emerges on the local scale in the real spectra. The statistical 
behavior of eigenvalues separated by distance S — sA measured in units of the mean eigenvalue 
spacing A is dictated by the global matrix symmetries (e.g. if they are complex Hermitian or real 
symmetric Q), being the same for all random matrix ensembles within a fixed symmetry class. 
All ensemble specific information is encoded in A. On different levels of rigor, this universality 
was established for "invariant" ensembles (i.e. matrices with invariant probabiltity distributions) 
p6|-p8[ and for matrices with i.i.d. entries, including sparse matrices [ p4| , ^ . Similar universality 
holds on a larger scale 5 ^ A p0| , |6]| and in the vicinity of the spectrum edges 

It turns out, that it is the local scale universality that is mostly relevant for real physical systems 
[Ql . Namely, statistics of highly excited bound states of closed quantum chaotic systems of quite 
different microscopic nature turn out to be independent of the microscopic details when sampled 
on the energy intervals large in comparison with the mean level separation, but smaller than the 
energy scale related by the Heisenberg uncertainty principle to the relaxation time necessary for a 
classically chaotic system to reach equilibrium in phase space Moreover, these statistics turn 
out to identical to those of large random matrices on the local scale, with different symmetry classes 
corresponding to presence or absence of time-reversal symmetry. 

One of the aims of the present paper is to demonstrate that complex spectra of weakly non- 
Hermitian random matrices possess a universality property which is as robust as the above men- 
tioned local scale universality in the real spectra of Hermitian matrices. Weakly non-Hermitian 
matrices appear naturally when one uses the Heidelberg approach to describe few-channel chaotic 
scattering [p^ . When the number M of open channels is small in comparison with the number N 
of the relevant resonances, the majority of the S-matrix poles (resonances) are situated close to the 
real axis. This is well captured within the Heidelberg approach. With a proper normalization of H 
and W, the imaginary part of typical eigenvalues of the effective Hamiltonian Tieff is of the order 
of the mean separation between neighboring eigenvalues along the real axis. This latter property 
is a characteristic feature of the regime of weak non-Hermiticity. 

Motivated by this example we introduced in |^5| another ensemble of weakly non-Hermitian ran- 
dom matrices. This ensemble consists of almost-Hermitian matrices which interpolate between the 
Gaussian ensemble of Hermitian matrices (GUE) and the Gaussian ensemble of complex matrices 
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studied by Ginibrc. It turned out that the eigenvalue distribution for ahiiost-Hermitian random 
matrices is described by a formula containing as two opposite limit cases both the Wigner 
semicircular distribution of real eigenvalues and the uniform distribution of complex eigenvalues 
obtained by Ginibre. Further studies of almost-Hermitian random matrices showed that ac- 
tually all their local scale eigenvalues statistics describe crossover between those of the GUE and 
Ginibre ensembles. Later on Efetov, in his studies of directed localization ||2^, discovered that 
weakly non-Hermitian matrices are relevant to the problem of motion of flux lines in superconduc- 
tors with columnar defects. Efetov's matrices are real almost-symmetric. They interpolate between 
Gaussian ensemble of real symmetric matrices (GOE) and the Gaussian ensemble of real asymmet- 
ric matrices. This development clearly shows that, apart from being a rich and largely unexplored 
mathematical object, weakly non-Hermitian random matrices enjoy direct physical applications 
and deserve a detailed study. 

The present paper consists of two parts. In the first part we study a three parameter family of 
random matrix ensembles which contains the above mentioned ensembles of almost-Hermitian and 
almost-symmetric matrices. Our random matrices are of the form 

H = {Si+ iuAi) + iv{S2 + iwA2), 

where the four matrices on the right-hand side are mutually independent, with 5i.2 being real 
symmetric and Ai_2 being real skew-symmetric. By choosing matrix distributions and varying 
the parameter values one obtains different ensembles of non-Hermitian matrices. We use that 
normalization of matrix elements which ensures that 

lim i-(TrS'^)== lim ^{TtA.AJ) = 1, j = 1,2 

N being the matrix dimension. The parameters v and v are scaled with matrix dimension: 

a 4> 
V = — u = — —, 
2//V 2Vn 

and a, (j), and w are assumed to be of the order of unity in the limit — > cx). The above scaling of 
V provides access to the regime of weak non-Hermiticity, while scaling u we describe the crossover 
between the GOE and GUE types of behavior of eigenvalues of the Hermitian part of H. A simple 
argument [^5| based on the perturbation theory shows that for our random matrices the eigenvalue 
deviations from the real axis are of the order of when N is large, i.e. it is of the same order as 
typical separation between real eigenvalues of the Hermitian Si + iuAi. Hence, in order to obtain 
a nontrivial eigenvalue distribution in the limit N oo one has to magnify the imaginary part 
scaling it with the matrix dimension. 

Our study of the scaled eigenvalues of H is based on the supersymmetry technique. We ex- 
press the density of the scaled eigenvalues in the form of a correlation function of a certain zero- 
dimensional non-linear cr-model. The obtained correlation function is given by a supersymmetric 
integral which involves only the density of the limit eigenvalue distribution of the Hermitian part 
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of H and the parameters a, (j), w. In two particular cases this supersymmetric integral can be 
explicitly evaluated yielding the earlier obtained distributions of complex eigenvalues for almost- 
Hcrmitian and almost-symmetric matrices [ p7| . 

The supersymmetric cr-model was invented long ago by Efetov in the context of theory of disor- 
dered metals and the Anderson localization and since then have been successfully applied to diverse 
problems [ |6^j65| . Application of this technique to the calculation of the mean density of complex 
eigenvalues of non-Hermitian random matrices was done for the first time in our earlier works 
|p^ , ^j45| and further advanced by Efetov |2^] in the context of description of flux line motion in 
a disordered superconductor with columnar defects. 

A detailed account of our calculations is given for sparse matrices ]5^-|55| with i.i.d. entries, 
although our results are extended to "invariant" ensembles and conventional random matrices 
with i.i.d. entries. We assume that matrix entries of Sk and Ak are distributed on the real axis 
with the density 

V{x) = (l - f ) Six) + ^h{x), (2) 

where h[x) is arbitrary symmetric density function, h{x) = h{—x), having no delta function sin- 
gularity a,t X — and satisfying the condition J x'^h{x)dx < oo. We also assume that the mean 
number of nonzero entries p exceeds some threshold value: p > pi, see ]59| . 

We want to stress that Eq. (H) describes the most general class ofrandom matrices whose entries 
are i.i.d. variables with finite second moment In particular, in the first part of our paper we 
do not assume the matrix entries to be Gaussian. 

We believe that here the power of the supersymmetry method is the most evident and we are 
not aware of any other analytical technique allowing to treat this general case non-perturbatively. 

Although giving an important insight into the problem, the supersymmetry non-linear cr— model 
technique suffers from at least two deficiencies. The most essential one is that the present state 
of art in the application of the supersymmetry technique gives little hope of access to quantities 
describing correlations between different eigenvalues in the complex plane due to insurmountable 
technical difficulties. At the same time, conventional theory of random Hermitian matrices sug- 
gests that these universal correlations are the most interesting features. The second drawback is 
conceptual: the supersymmetry technique itself is not a rigorous mathematical tool at the moment 
and should be considered as a heuristic one from the point of view of a mathematician. 

In the second part of the present paper we develop the rigorous mathematical theory of weakly 
non-Hermitian random matrices of a particular type: almost-Hermitian Gaussian. Our consider- 
ation is based on the method of orthogonal polynomials. Such a method is free from the above 
mentioned problem and allows us to study correlation properties of complex spectra to the same 
degree as is typical for earlier studied classes of random matrices. The results were reported earlier 
in a form of Letter-style communication Unfortunately, the paper contains a number 

of misleading misprints. For this reason we indicate those misprints in the present text by using 
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footnotes. 



II. REGIME OF WEAK NON-HERMITICITY: UNIVERSAL DENSITY OF COMPLEX 

EIGENVALUES 

To begin with, any N x N matrix J can be decomposed into a sum of its Hermitian and skew- 
Hermitian parts: J = Hi +ii?2, where Hi — {J + J^)/2 and H2 — {J — ,V)/2i. Following this, we 
consider an ensemble of random N x N complex matrices J = Hi + ivH2 where Hp] p — 1,2 are 
both Hermitian: Hj^ = Hp. The parameter v is used to control the degree of non-Hermiticity. 

In turn, complex Hermitian matrices Hp can always be represented as Hi — Si + iuAi and 
H2 = S2 + iwA2, where Sp — Sp is a real symmetric matrix, and Ap — —A^ is a real antisymmetric 
one. From this point of view the parameters u, w control the degree of being non-symmetric. 

Throughout the paper we consider the matrices Si, 82, Ai, A2 to be mutually statistically inde- 
pendent, with i.i.d. entries normalized in such a way that: 

lim ^TrSl = lim ^TrApA^ = 1 (3) 

As is well-known [Q, this normalisation ensures, that for any value of the parameter u 7^ , 
such that u ~ 0(1) when iV — s- 00, statistics of real eigenvalues of the Hermitian matrix of the 
form H = S + iuA is identical (up to a trivial rescaling) to that of u = 1, the latter case known 
as the Gaussian Unitary Ensemble (GUE). On the other hand, for u = real eigenvalues of real 
symmetric matrix S follow another pattern of the so-called Gaussian Orthogonal Ensemble (GOE). 

The non-trivial crossover between GUE and GOE types of statistical behaviour happens on a 
scale u (X ]6^ . This scaling can be easily understood by purely perturbative arguments 

jSTf . Namely, for u cx \/N^/'^ the typical shift 5\ of eigenvalues of the symmetric matrix S due to 
antisymmetric perturbation iuA is of the same order as the mean spacing A between unperturbed 
eigenvalues : (5A~A~l/iV. 

Similar perturbative arguments show ]45| ] , that the most interesting behaviour of complex eigen- 
values of non- Hermitian matrices should be expected for the parameter v being scaled in a similar 
way: v cx 

It is just the regime when the imaginary parts Im of a typical eigenvalue Zk 
due to non-Hermitian perturbation is of the same order as the mean spacing A between unper- 
turbed real eigenvalues : IvaZk ~ A ~ Under these conditions a non-Hermitian matrix J still 
"remembers" the statistics of its Hermitian part Hi. As will be clear afterwards, the parameter w 
should be kept of the order of unity in order to influence the statistics of the complex eigenvalues. 

It is just the regime of weak non-Hermiticity which we are interested in. Correspondingly, we 
scale the parameters as 0: 



^In the Letter Ml[ there is a misprint in the definition of the parameter a. 



2VN 

and consider a,(f),w fixed of the order 0(1) when N oo. 
One can recover the spectral density 

N N 



(4) 



(5) 



k=l 



k=l 



of complex eigenvalues Zk = Xfc + iYk, k = 1,2, .... N from the generating function (cf. ||l^ , po|| ) 



Det[{Z~J){Z-jy + K'^] 
Dei[(Zfc- J)(Zb- J)t+K2] 



(6) 



as 



lini ^ „ lim -ttt^Z. 



■K K^o dZ* Zt^z dZf, 

To facilitate the ensemble averaging we first represent the ratio of the two determinants in Eq. (|^) 
as the Gaussian integral 

N 

I l[[d^^] exp{/:i($) + £2($)} (7) 

over 8-component supervectors 



i?»(±) 



x.(±) 



with components ?'i(+), 7'i(— ); « = 1,2,..., N being complex commuting variables and 
Xi(+), Xi{^) forming the corresponding Grassmannian parts of the supervectors 4'i(±). The terms 
in the exponent of Eq. (^ are of the following form: 



(8) 
(9) 



Here /2 = diag(l,l), /i = diag(/2,/2), A = diag(/4, -14), A^ = diag(/2, 62, -/2, 62), A/ = A - A^, 
S = + Ey, and 
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and the matrices "SrM, "^tJ, S-r, A^ are obtained from the corresponding matrices without subindex 
T by replacing all I2 blocks with the matrices f = diag(l, —1). 



We also use ^J(±) = (x*(±); — Xi(±)). When writing £2{^) in Eq.(||) we have used the fact that 
diagonal matrix elements S2.U for i — 1, .., N give total contribution of the order of 0{1/N) with 
respect to the total contribution of the off-diagonal ones and can be safely disregarded. 

Now we should perform the ensemble averaging of the generating function. We find it to be 
convenient to average first over the distribution of matrix elements of the real symmetric matrix 
Si. 

These elements are assumed to be distributed according to Eq.(||). Before presenting the deriva- 
tion for our case, let us remind the general strategy. The procedure consists of three steps. First 
step is the averaging of the generation function over the disorder. It can be done trivially due to 
statistical independence of the matrix elements in view of the integrand being a product of expo- 
nents, each depending only on the particular matrix clement Hij. This averaging performed, the 
integrand ceases to be the simple Gaussian and thus the integration over the supervectors can not 
be performed any longer without further tricks. When matrix elements are Gaussian-distributed, 
this difficulty is circumvented in a standard way by exploiting the so-called Hubbard- Stratonovich 
transformation. That transformation amounts to making the integrand to be Gaussian with re- 
spect to components of the supervector by introducing new auxilliary integrations. After that the 
integral over supervectors can be performed exactly, and remaining auxilliary degrees of freedom 
are integrated out in the saddle-point approximation justified by large parameter N. 

As is shown in the paper |]5^ , there exists an analogue of the Hubbard-Stratonovich transforma- 
tion allowing to perform the steps above also for the case of arbitrary non-Gaussian distribution. 
The main difference with the Gaussian case is that the auxilliary integration has to be chosen in 
a form of a functional integral . 

Our presentation follow the procedure suggested in , and presented also in some detail in |55| 



Exploiting the large parameter N ^ 1 one can write: 



(exp £!($)) |jv>i« exp 



P 

2N 



; hp{z) = / dsh{s)e 



(10) 



In order to proceed further we employ the functional Hubbard-Stratonovich transformation 
introduced in IM: 



exp 



P 

2N 



2J 



(11) 



Vg exp 



pN 



N 



dedeg{e)c{e, 6)5(6) -f p ^ g{<^>,) 



■^Similarly to the paper jwj we first disregard necessity for the compactification and use the matrix A 
rather than two different matrices L and A, see discussion in 
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where the kernel C(8, 9) is determined by the relation: 

J dec(e, e)/iF(eA$) = s{q, $) (12) 

with the right-hand side of the eq.(p^ being the (5— function in the space of supervectors. 

Substituting eg. (pi]) into averaged eq.(0) and changing the order of integrations over and 
Vg one obtains the averaged generating function in the form: 

(Z) ^Jvg exp [-N£{g) + S£{g)] (13) 



where 



^9) = 1/ mm ff(<i>)c($,$)5(i>) -In J m] e^(*) 

IUllm^] expE.^(<i>.)] 



with T{<P) = §X$TA$ + p.g($), 
1 " 

and f ^Ki~i{X ~ Xb)Ab - nSr,b - 

We are interested in evaluating the functional integral over Vg in the limit N oo and k —^ 
0; X Xb- Moreover, we expect eigenvalues of weakly non-Hermitian matrices to have imaginary 
parts Y to be of the order of Remembering also the chosen scaling (^, we conclude that 

the argument of the logarithm in Eq.([l4|) is close to unity and the term SC{^) in Eq.(p^ should 
be treated as a small perturbation to the first one. Then the functional integral of the type 
J 'Dg{...) exp —NC{g) can be evaluated by the saddle-point method. Variating the "action" C{g) 
and using the relation eq.(|l^) one obtains the following saddle point equation S£{g)/Sg = for 
the function .g($): 

/[d9]fef(<l>tL9)exp.F(9) 
^('^^ = /[d9]exp.F(e) 



A quite detailed investigation of the properties of this equation was performed in 1 54 55[| . Below 
we give a summary of the main features of the eq.(^5|) following from such an analysis. 

First, the solution .g($) to this equation can be sought for in a form of a function gs{^) = go{x, y) 
of two superinvariants: x = <i>^^<i> and y = <i>^^A$. 



As the result, the denominator in eq.(15) is equal to 1 due to the identity J[d^]F{x, y) = F(Q, 0) 
which is a particular case of the so-called Parisi-Sourlas-Efetov-Wegner (PSEW) theorem, see e.g. 
|l65| and references therein. However, the form of the function go{x,y) is essentially different for 
the number of nonzero elements p per matrix column exceeding the threshold value p ~ pi and for 
p < pi [5^j59[ |. Namely, for p > pi the function gQ{x, y) is an analytic function of both arguments x 
and y, whereas for p < pi such a function is dependent only on the second argument y = <i>^A$. 
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At the same time, the saddle-point equation eq.(|15|) is always invariant w.r.t. any transformation 
.g($) g{T^) with supermatrices T satisfying the condition T'^AT = A. 

Combining aU these facts together one finds, that for p > pi a. saddle-point solution .gs($) gives 
rise to the whole continuous manifold of saddle-point solutions of the form: gri^) = gs{T^) = 
go{^''TT^,^''A^), so that all the manifold gives a nonvanishing contribution to the functional 
integral eq-dl^). It is the existence of the saddle-point manifold that is actually responsible for the 
universal random-matrix correlations . 

We see, that the saddle-point manifold is parametrized by the supermatrices T. It turns out, 
however, that one has to "compactify" the manifold of T matrices with respect to the "fermion- 
fermion" block in order to ensure convergence of the integrals over the saddle-point manifolds [ p2[ . 
The resulting " compactified" matrices T form a graded Lie group U0SP(2, 2/4). Properties of 
such matrices can be found in [ p2[ together with the integration measure dpL(T). 

^From now on we are going to consider only the case p > pi. The program of calculation is as 
follows: (i) To find the expression for the term SC{g) on the saddle-point manifold g — gri^) in 
the limit — s- oo and (ii) to calculate the integral over the saddle-point manifold exactly. 

Expanding the expression in Eq. (|l4|) to the first non- vanishing order in /, 82-, Ai, A2, introducing 
the notation J^t{^) — ^X^^A^ + pgxi^) and using the relations 



jme^pTrm = 1; JULiid^k] ( $15$, ) |.<, exp ^T(<i>/c) - 



which hold for arbitrary 8x8 supermatrices S, (7, one finds that 

5C{gT) = \{^^f^)T + SCn + SCjn, 

where 



(16) 



(17) 



U 

~2 



2 2 



iuvw 



i<j 



A2 



${A^$2 



${S$2 



and we used the notations 



[d$i][d$2] ($Ii?$2) (*IC'$2) exp[J^T(*i) +^t($2 
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It is clear that with the chosen normaHzation [see Eqs. - (Q)] we have 



i<j i<j 

when N oo. On the other hand, it is easy to see that: 



i<j 



2 u;2^2 
ij 8 



(18) 



(19) 



because of the statistical independence of 52,^1,^2 and the chosen normalization of matrix ele- 
ments. Therefore, the part SCj^i can be safely neglected in the limit of large N. 

To proceed further it is convenient to introduce the 8x8 supermatrix W with elements 



(20) 



Exploiting the saddle-point equation for the function gri^) = go ^'I'^T^T^, $^A$^ one can show 
(details can be found in Eqs. (67-70)) that the supermatrix W whose elements are Wafs can 
be written as: 



-4 



goyA + igoxAQ 



(21) 



where B is the second moment of the distribution h{s): B = J h{s)s^ds and gox = 
dgo/dx\x=y=o', goy ~ dgQ/dy\x=y=o- In the expression above we introduced a new supermatrix: 
Q = -if'^kf. 

Using the definition of the matrix W, one can rewrite the part SCr as follows (cf. [Q, eqs. 
(71)-(73)): 



Str WY.rWJ:r - (f>^ Str WArWAr + w^(f>^ Str WY.W^ 



(22) 



Now one can use Eq.(20) together with the properties: StrQ = Str A = Str / = 0; ~ —I to 
show that: 



9ox 
4B2 



Str Q&rQc^T — 0^ Str QT2QT2 + w'^cj)'^ Str Q&Qa 



















I 















where the 8x8 supermatrices entering these expressions are as follows: 

T2 = diag{f3,f3}; ct^ = 

and fa is 4 X 4 diagonal, fa = diag{'r, f }. 
In the same way one finds: 



' T B 

2igo: 



B 



K Str QA - i{X - Xb) Str KbQ - Is Str ct^^'Q - Y Str a^^^Q 
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where 

i-B = diag{/2, 62,72,02}; 



^3 



and f^'^ are 4x4 diagonal supermatrices: t^^'' = diag{f,02} and fg^^ = diag{02,T}. 

At last, we use the relation between gQ{x,y) and the mean eigenvalue density for a sparse 



symmetric matrix Si at the point X on the real axis derived in |54|: 



:y(X)^^{p{X))=~^go. (23) 

Substituting expressions for SCr and to the generating function Z represented as an 

integral over the saddle-point manifold parametrized by the supermatrices T (or, equivalently, by 
the supermatrices Q — —iT^^KT) and performing the proper limits we finally obtain: 

(p(X, Y)) = j d^i{Q) Str (af )q) Str (a,g) cxp -S{Q) (24) 

5(Q) = -^yStr(<7,Q) -^Str(<7,Q) + _ Str (tsQ) - ^ Str (^q) 

where we introduced the scaled imaginary parts y — 7n'{X)NY and used the notations: = 
{Triy{X)af , = {Trv{X)(l)f , = {7rv{X)awf 

The expression ( p^ is just the universal a~ model representation of the mean density of complex 
eigenvalues in the regime of weak non-Hermiticity we were looking for. The universality is clearly 
manifest: all the particular details about the ensembles entered only in the form of mean density 
of real eigenvalues i^iX). The density of complex eigenvalues turns out to be dependent on three 
parameters: a,b and c, controlling the degree of non-Hermiticity (a), and symmetry properties of 
the Hermitian part (b) and non-Hcrmitian part (c). 

The following comment is appropriate here. The derivation above was done for ensembles with 
i.i.d. entries. However, one can satisfy oneself that the same expression would result if one start 
instead from any "rotationally invariant" ensemble of real symmetric matrices Si. To do so one can 
employ the procedure invented by Hackenbroich and Weidenmiiller allowing one to map the 
correlation functions of the invariant ensembles (plus perturbations) to that of Efetov's cr— model. 

Still, in order to get an explicit expression for the density of complex eigenvalues one has to 
evaluate the integral over the set of supermatrices Q. In general, it is an elaborate task due to 
complexity of that manifold. 

At the present moment such an evaluation was successfully performed for two important cases: 
those of almost-Hermitian matrices and real almost-symmetric matrices. The first case ( which is 
technically the simplest one) corresponds to — + cx3, that is — > 00. Under this condition only that 
part of the matrix Q which commutes with T2 provides a nonvanishing contribution. As the result, 
Str (^a-Qj — Str (^arOj so that second and fourth term in Eq.(p^) can be combined together. 
Evaluating the resulting integral, and introducing the notation = + one finds Ea]: 
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2 1 

Px{y) - Y^^exp (^-^) J dtcosh{2ty)exp{-&Hy2), (25) 



where px{y) is the density of the scaled imaginary parts y for those eigenvalues, whose real parts 
are situated around the point X of the spectrum. It is related to the two-dimensional density as 
px{y) = p{X,Y)/7T{Ni.{X))\ 

It is easy to see, that when a is large one can effectively put the upper boundary of integration 
in Eq. ( |25| ) to be infinity due to the Gaussian cut-off of the integrand. This immediately results in 
the uniform density px{y) — (5^)^^ inside the interval \y\ < and zero otherwise. Translating 
this result to the two-dimensional density of the original variables X, Y, we get: 

1 .-^7? ?T ^i\Y\<2Tn^{Xy{l + w^) 

piX,Y) = { Anv^l + w^) (26) 

otherwise 

This result is a natural generalisation of the so-called "elliptic law" known for strongly non- 
Hermitian random matrices ||36| , [r^ . Indeed, the curve encircling the domain of the uniform eigen- 
value density is an ellipse: 2v^{i+w'^) ~^ ^ = 1 as long as the mean eigenvalue density of the 
Hcrmitian counterpart is given by the semicircular law, Eq.(0) (with the parameter J = 1). The 
semicircular density is known to be shared by ensembles with i.i.d. entries, provided the mean 
number p of non-zero elements per row grows with the matrix size as p cx N"; a > 0, see [ ^4[ . 
In the general case of sparse or "rotationally invariant" ensembles the function I'iX) might be 
quite different from the semicircular law. Under these conditions Eq. ( p6| ) still provides us with the 
corresponding density of complex eigenvalues. 

The second nontrivial case for which the result is known explicitly is due to Efetov . It is the 
limit of slightly asymmetric real matrices corresponding in the present notations to: (/> — > 0; w — > oo 
in such a way that the product (l)w = c is kept fixed. The density of complex eigenvalues turns out 
to be given by: 

1 

Px{y) = S{y) J dtcxp{^£hy2) (27) 



+ 2y^M J duexp(^-?^^^ J dttsmh{2t\y\)exp{-d'f/2), 

1 

The first term in this expression shows that everywhere in the regime of "weak asymmetry" 
c < oo a finite fraction of eigenvalues remains on the real axis. 

Such a behaviour is qualitatively different from that typical for the case of "weak non- 
Hermiticity" a < oo, where eigenvalues acquire a nonzero imaginary part with probability one. 

In the limit c >> 1 the portion of real eigenvalues behaves like c^^. Remembering the normali- 
sation of the parameter v, Eq.(|^), it is easy to see that for the case of v — 0(1) the number of real 
eigenvalues should scale as \/iV. Indeed, as was first noticed by Sommers et al. the number 
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of real eigenvalues of strongly asymmetric real matrices is proportional to a/TV . This and the fact 
that the mean density of real eigenvalues is constant was later proved by Edelman et al. |^8|. 

III. GAUSSIAN ALMOST-HERMITIAN MATRICES: FROM WIGNER-DYSON TO 

GINIBRE EIGENVALUE STATISTICS 

In the previous section we obtained the eigenvalue distribution in the regime of weak non- 
Hermiticity for the random matrices of the form J — Hi + ivH2 , with Hi and H2 being mutually 
independent Hermitian random matrices with i.i.d. matrix entries. The obtained eigenvalue distri- 
bution appeared to be universal, i.e. independent of the probability distribution of the Hermitian 
matrices Hi and i?2- 

In the present section we reexamine a particular case of J when both Hi and H2 are taken to 
be Gaussian. In this special case not only the mean eigenvalue density but also the eigenvalue 
correlation functions can be obtained and studied in great detail. 

The ensemble of random matrices that will be considered in this section is specified by the 
probability measure dfj,{J) — V{J)dJ, 

(-^)" exp[-^Tr(iit_KeP)] (28) 
on the set M. of complex N x N matrices with the matrix volume element 

N 

dJ = Y[ d^Jjk, d^Jjk = dKe JjkdlmJjk- 

If the Hermitian Hi and H2 are taken independently from the GUE, the probability distribution 
of J = Hi + ivH2 is described by the above-given measure dfi{J) with 



1 + ^2 

provided that Hi and H2 are normalized to satisfy {Tr H^) = N{1 + r)/2, p ^ 1,2. 

The parameter t, < r < 1, controls the magnitude of the correlation between Jjk and Jkj'- 
(JjkJkj) = t/N, hence the degree of non-Hermiticity. All Jjk have zero mean and variance 
{\Jjk\'^) — 1/-^ and only J^fe and J^j are pairwise correlated. If t = all Jj^ are mutually 
independent and we have maximum non-Hermiticity. When r approaches unity, Jj^ and J^^ are 
related via Jju — J^j and we arc back to the ensemble of Hermitian matrices. 

Our first goal is to obtain the density of the joint distribution of eigenvalues in the random matrix 
ensemble specified by Eq. (p8|). First of all, one can disregard the matrices whose characteristic 
polynomial has multiple roots. For, the set of such matrices forms a surface in Al, hence has zero 
volume. Every matrix off this surface has iV distinct eigenvalues and we label them Zi, . . . , Zjv 
ordering them in such a way that 

\Zi\<\Z2\<...<\Zn\ (29) 
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and if \Zj\ = \Zj^i\ for some j then argZj < argZj^i. Given J, one can always find a unitary 
matrix U and a triangular T (Tjk = if j > /c) such that 



(30) 



and Tjj = Zj for every j [ pSj . The choice of U and T is not unique. For, multiplying U to the 
right by a unitary diagonal matrix $ one can also write J = VSV~^, where V = is unitary, 
S is triangular, and again Sjj — Zj for every j. It is natural, therefore, to impose a restriction on 
U requiring, for instance, the first non-zero element in each column of U to be real positive. Then 
the correspondence ( ^0|) between J and {U,T) is one-to-one. 

The idea of using the decomposition (js^) ( which is often called the Schur decomposition) for 
derivation of the joint distribution of eigenvalues goes back to Dyson |l69| and we simply follow 
his argument. To obtain the density of the joint distribution one integrates ( p8[ ) over the set of 
matrices whose eigenvalues are Zi, . . . , Zj^. To perform the integration, one changes the variables 
from J to {U, T) and integrates over U and the off-diagonal elements of T. The Jacobian of the 
transformation J —> ([/, T) depends only on the eigenvalues and is given by the squared modulus 
of the Vandermonde determinant of the {Zj} |6^. Since IJ is unitary. 



Tr(ij1' +TReP) = Tr{ff''+TRef^). 
Therefore, the integral over U yields 

Vol[;7(iV)] 7r^(^-i)/2 



(31) 



(27r) 



N 



N-1 



n 



where Yo\[U{N)] is the volume of the unitary group U{N). Since T is triangular, the integration 
over the off-diagonal entries of T reduces, in view of (|3l|), to the Gaussian integral 



Y[ I d^Tjk exp 



N 



1 — T' 



(|r,,|2_rRe7ffe) 





'nil -T^y 




N 



l<j<k<N ' 

Collecting the constants one obtains the desired density. Obviously, it is symmetric in the eigen- 
values {Zj}. Therefore, the above restriction of Eq. ( p^ ) on the eigenvalues can be removed by 
reducing the obtained density in A^! times. Thus finally, the density of the joint distribution of 
(unlabelled) eigenvalues in the random matrix ensemble specified by Eq. ( |2^ ) is given by 

N 

Vn{Z,, ...,Z„)= IT w\Z,) W \Z, Z,\\ (32) 



^^(^+1)72 

^^ ii...iv'(i-r^w2 n ^'(^^o n 1^^- - ^'=1'' 



where 



w^{Z) = exp<^ - 



N 



\Z\ 



z* 



(33) 



The form of the distribution Eq. ( |3^ ) allows one to employ the powerful method of orthogonal 
polynomials H. Let Hn{z) denotes nth Hermite polynomial. 
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Hniz) = -^^^exp J dt t'"- exp f-~L.zfizt]. 



(34) 



These Hermite polynomials are orthogonal on the real axis with the weight function cxp(— a;^/2) 
and are determined by the following generating function 

+2^ 



exp ( zi - ^) = ^ Hn{z)- 



It is convenient to rescale Hermite polynomials in the following way: 

^ ' 0FV^(1 - r2)i/4 " I V T i 



(35) 



(36) 



The main reason for doing that rescaling is that these new polynomials n = 0, 1, 2, . . ., are 

orthogonal in the complex plane Z = X + iY with the weight function w^(Z) of Eq. (^3|): 



d'Z p,,{Z)p^{Z*)w\Z) = 5, 



(37) 



(Recall that d^Z = dXdY.) We borrowed this observation which is crucial for our analysis from 
the paper (see also the related paper |^). A quick check of the orthogonality relations is 
possible with the help of the generating function (^). 

With these orthogonal polynomials in hand, the standard machinery of the method of orthogonal 
polynomials |^ yields the n-eigenvalue correlation functions 

iV! 



i?„(Zi, .... Z„) — — — — — - / d^Zn+i---d^ZN Vn{Z} 
[N - ny. J 



in the form 



Rn{Zi, Z„) = det [K]y{Zj, Zk)]2k=i ^ 
where the kernel Kn{Zi, Z2) is given by 

K^iZi, Z2) - W{Z,)W{Z^) Pn{Zl)Pn{Z^) 



(38) 



(39) 



(40) 



In particular, define the density of eigenvalues as in Eq.(||), so that the number of eigenvalues in 
domain A of the complex plane is given by the integral 



n{A) = J d^Z p{Z). 

A 



(41) 



Notice that the averaged density of eigenvalues {p{Z)) is simply Ri{Z). From Eqs. (p9[), (|4C|), and 
( ^ ) one infers that 

2 



Ri{z) 



N 



exp 



N 



1 -t2 



\z\^--{z^ + z- 



,2^ I I ^ 



y - 



(42) 
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This exact result is valid for every finite N. The rest of this section is devoted to sampling 
of statistical information that can be obtained from Eqs. (|39| ) - ( ^ ) for large matrix dimensions 
iV 3> 1. First we briefly examine the regime of strong non-Hermiticity when the real and imaginary 
parts of a typical eigenvalue are of the same order of magnitude when — > oo. This regime is 
realized when < limTv^oo t < 1 (recall that for r = 1 our matrices are Hermitian). We will show 
that in this case r-dependence of the eigenvalue correlations on the local scale becomes essentially 
trivial and the correlations become identical to those found by Ginibre in the case of maximum 
non-Hermiticity {t — 0). 

Then we will examine the regime of weak non-Hermiticity when the imaginary part of typical 
eigenvalue is of the order of the mean separation between the nearest eigenvalues along the real 
axis. This regime is realized when 

T=l-|^, a>0. (43) 

We will show that by varying the parameter a one can describe the crossover from the Wigner- 
Dyson eigenvalue statistic typical for Hermitian random matrices to the Ginibre eigenvalue statistic 
typical for non-Hermitian random matrices. 

To begin with the regime of strong non-Hermiticity we first recall that in this regime the eigen- 
values in bulk are confined to an ellipse in the complex plane and they are distributed there with 
constant density (cf. Eq.(p6|)) : 

r 1 . 

lin, i-i?i(Z)=i ^{l-T^y " (l + r)2 + (l-r)2 - 
0, otherwise. 



This fact can be inferred from Eq. ([421). Inside the ellipse every domain of linear dimension of the 
order of contains, on average, a finite number of eigenvalues. Thus the eigenvalue statistics 

on the local scale are determined by the correlation functions 



i?„(zi, . . . , z„) = iV-"i?„(ViVZi, . . . , y/NZn 



of the rescaled eigenvalues z — \lNZ. This rescaling is effectively equivalent to the particular 
normalization of the distribution (^8|) which yields (Tr JJ^^ — iV^, the normalization used in [Q. 
One can easily evaluate the rescaled correlation functions Rn exploiting Mehler's formula^ [ f72[ : 

|'S«..(^)«..(^) =^»p{t^[-:-5(^?-")]}- 

Indeed, denote K]ss{z\^Z2) = N^'^Ki^{\rN zx^-JN Z2)- Then, by Mehler's formula 



^Mehler's formula can be derived by using the integral representation (^) for Hermite polynomials in 
the Lh.s. of Eq. (||). 
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lim Kn{zi,Z2) = 



1 



N^oo - '■ ' 7]-(l — r^) 
and in view of the relationship 



exp 



2^122* -|Z1|2-|Z2|2 



one obtains that 



N 



hm Rn{zi, . . . ,z„) = 



2(1 -r2) 
i?„(zi, ...,z„) = det kN{zj,Zk) 



exp 



T{zt 



zl + zl-zf) 



4(1 -r2) 



7r(l-T2 



dot 



exp 



l-r2 



In particular 



lim Ri{z) = — — 

N^oc 7r(l - r^) 



N 



lim i?2(2i, Z2) = 



1 



n 2 



7r(l -r2) 



exp 



1-t2 



(45) 

(46) 
(47) 



After the natural additional rescaling z/^1 — t"^ z Eqs. ( p5| ) - ( pTf ) become identical to those 
found by Ginibre |^6| for the case r = 0. 

Now we move on to the regime of weak non-Hermiticity [see Eq. (|4^)]. We will show that in this 
regime new non-trivial correlations occur on the scale: \a\Z\2 = 0{l/N), ReZi— RcZ2 = 0{1/N). 

To find the density of complex eigenvalues and describe their correlations, let us define new 
variables yi, y2, oj: 

Zi = A H h ? — , Z2 = X hi — . (48) 

Our first goal is to evaluate the kernel Kn{Zi, Z2) in the limit 

N ^ oo,N{1-t) ^ X, y 1,2 are fixed. (49) 

Using the integral representation for Hermite polynomials, Eq. (^), we can rewrite Kn{Zi, Z2) in 
the form 



Kn{Zr,Z2) 



N 



■n\J\ — T' 



: exp 



NiZl^^Zl'') A^(|Zi|2 + |Z2*|2) Nt{ZI^ Zl + Zl^ ^ Z. 



2t 



^^2^ /ds(is)"exp 



— 00 —00 



2(1 -r2) 



f2 + s2 



4(1 -r2 



{itZi — 13^2) 



Using new variables (Eq) in the equation above and making the substition 



u — (t + s)\ — and v — (t — s)\ — 
V r V r 



in the integrals, we obtain 

iV2 



KNiZuZ2) = 



: exp 



V1+V2 , iX{yi-V2) , iuj{vi+V2) 



rVl - t2 [4iVT(l + T) 2A^t(1-t) r 27Vr 

iXvN NX^ 1 i;(yi - ^2) 



1 
2^ 

1 

2^ 



—00 



dv exp 
du exp 



Nv^ / 1 



7Vm2 



T t(1+t) 

mti; u{yi + ^2) 



2r 



2t 



2t 



6 



TV 



(u2 - i;2) 
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were we have introduced the notation 



N-l 



QNix) 



n=l 



Now we are in a position to evaluate Kn{Zi, Z2) in the regime defined by Eq.(|4S|). Indeed, in 
this regime, to the leading order, 

iXvN NX"^ 1 _ N{v-iX) 









r 4 


-I) 


f Nu^ 




-0 


I 4 





T t{1+t) 

iuuj _ u[yi + 112) 
2t 



2t 



iP'o? iuuj u{yi — 1/2) 



i From these relations one obtains that 



1 

Kn{Zi,Z2) - exp 

ZTT a 



du 

7^ 



exp 



u^a^ iuuj u{yi + y2) 



N 



dv exp 



2 2 
iV(w - v{yi - 2/2) 



e 



TV 



N 



Taking into the account that 



lim QNiNx) 

TV— 00 



1, < a; < 1 
0, a; > 1 



(50) 



and evaluating the integral over v by the saddle point method we finally obtain that in the regime 
Eq. (ES), to the leading order. 



Kn I X 



w , iy2 



2N N ' 2N N 
where y = {yi + ?/2)/2 and 

TTV(X) 

9a{.y) j 



7ra 



du 



■ exp 



2/? + 2/1 , iX{yT,-y2) 



exp 



2 7 
a u 



2uy 



.9" ( y - y ) > (51) 



(52) 



with v{X) = 2^\/4 — X'^ standing for the Wigner semicircular density of real eigenvalues of the 
Hermitian part Hi of the matrices J. 

Equation ( [5l| ) constitutes the most important result in this section. The kernel -ftTjv given by Eq. 
( pT| ) determines all the properties of complex eigenvalues in the regime of weak non-Hcrmiticity. 
For instance, the mean value of the density p{Z) ~ — Zi) of complex eigenvalues 

Z = X + iY is given by {p{Z)) = Kn{Z, Z). Putting yi = ya and w = in Eqs. ©-(H) we 
immediately recover the density Eq. ( p5| ) found by the supersymmetry approach |^ 



■^In the present section we normalized H2 in sucli a way tliat for weak non-Hermiticity regime we have 
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One of the most informative statistical measures of the spectral correlations is the 'connected' 
part of the two-point correlation function of eigenvalue densities: 



(p(Zi)p(Z2)), = <5(2)(Zi - Z2) - 3^2(^1, ^2), 



(53) 



In particular, it determines the variance S^(-D) = (n(Z?)^) — {n{D))'^ of the number n = 
J (PZp(Z) of complex eigenvalues in any domain D in the complex plane, see the Appendix for a 

D 

detailed exposition. 



Comparing with the definitions, Eqs. ( |3q ) and ( |39| ) we see that the cluster function 3^2(-^i, ^2) 
is expressed in terms of the kernel ifjv as 3^2(^1, ^2) = \Kn{Zi, ^2)!^- 

It is evident that in the limit of weak non-Hermiticity the kernel Kn in Eq.(|5l|) depends on 
X only via the semicircular density i^iX). Thus, it does not change with X on the local scale 
comparable with the mean spacing along the real axis A ^ 1/N. 

The cluster function 3^2 is given by the following explicit expression: 

2 



iV4 

3^2(^1, 2^2) = -^r^ exp 



2(2;? 



2/1) 



■!TU{X) 



du 



exp 



- u{yi + 2^2) + iuui 



(54) 



The parameter a = 'Kv{X)a controls the deviation from Hermiticity. When a ^ 1 the limits of 
integration in Eq.(|5^) can be effectively put to ±00 due to the Gaussian cutoff of the integrand. The 
corresponding Gaussian integration is trivially performed yielding in the original variables Zi, Zi 
the expression equivalent (up to a trivial rescaling) to that found by Ginibre |^^: y-iiZx^Z-i) = 
[N"^ /tto^Y exp{— A^^jZi — ^2^/0;^}. In the opposite case a ^ the cluster function tends to GUE 



form 3^2 (w,y 1,2/2) = -^5{yi)5{y2 



1 sin 7rz^(X)a; 



One can also define a renormalized cluster function: 

3^2(^1,^2) 



^"2(^1,^2) = 



Ri{Zi)Ri{Z2y 



Introduce the notation 



^2(1^, 2/1,2/2) = lim Y2 



X 



2N 



N ' 2N 



W2 



Then 



5^2(^^,2/1, 2/2) = 



dt exp ( — 



cosh [t{yi + y2 + iio)] 



dt exp ( — 



cosh(22;?/i) / dt exp — 



(55) 



cosh(2i2/2) 



limjv^cx) TrffI = N, whereas the normalization Eq.(^ gives limjv^oo Tr/}"! = A(l + w^). It is just because 
of this difference the parameter a entering Eq.(^5|) contains extra factor l + w"^ as compared to the present 
case. 
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where 

2/1,2 = yi,2T^v{X), Co = UJTTV{X) 

It has advantages of being non-singular in the Hcrmitian hmit a — > and coinciding with the 
usual GUE cluster function Cj^'^ sit? Co on the real axis: 2/1=2/2 = 0. We plotted this function for 
different values of the parameter a in Fig. |l|. 

The operation of calculating the Fourier transform of the cluster function over its arguments 
^1 2/1 1 2/2 amounts to simple Gaussian and exponential integrations. Performing them one finds the 
following expression for the spectral form-factor. 

00 oc 00 

b{qi,q2,k)^ j duj J dyi y (iy2 3^2(^1, ^2) exp[27ri(tjfc + yiqi + 2/292)] (56) 

— 00 ~OQ —OQ 

a^qj + ql + 2k^)] sin[TT^a^{qi + q2){l^{X) ~ \k\ ^ 



exp 



-%(X)-|fc|), 



7r2Q;2(gi + 92) 

where Zi and Z2 are given by Eq. (^8|), and diu) = 1 for u > and zero otherwise. 

We see, that everywhere in the regime of weak non-Hermiticity Q < a < 00 the formfactor shows 
a kink-like behaviour at |fc| = viX). This feature is inherited from the corresponding Hermitian 
counterpart-the Gaussian Unitary Ensemble. It reflects the oscillations of the cluster function with 
Lo which is a manifestation of the long-ranged order in eigenvalue positions along the real axis [Q. 
When non-Hermiticity increases the oscillations become more and more damped as is evident from 
the Fig. [|. 

As is well-known [||j4j, the knowledge of the formfactor allows one to determine the variance S2 
of a number of eigenvalues in any domain D of the complex plane. Small E2 is a signature of a 
tendency for levels to form a crystal-like structure with long correlations. In contrast, increase in 
the number variance signals about growing decorrelations of eigenvalues. 



For the sake of completeness we derive the corresponding relation in the Appendix, see Eq.(A3). 
In the general case this expression is not very transparent, however. For this reason we restrict 
ourselves to the simplest case, choosing the domain D to be the infinite strip of width (in 
units of mean spacing along the real axis A = (i/(0)iV)~^) oriented perpendicular to the real axis: 
< RcZ < L2,A; —00 < ImZ < 00. Such a choice means that we look only at real parts of 
complex eigenvalues irrespective of their imaginary parts. It is motivated, in particular, by the 
reasons of comparison with the GUE case, for which the function Y,2{Lx) behaves at large 
logarithmically: £2(^0;) oc InLj, [Q. 

After simple calculations (see Appendix) one finds ^ 



^In our earlier Letter the expression Eq.(^^ and formulae derived from it erroneously contained vra 
instead of a. 



S2(i.)=ix< 1-2 dk 1- 



sin^(7rfc) 
(7rfc)2 



■ exp 



ak 

L'l- 
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(57) 



First of all, it is evident that S2 grows systematically with increase in the degree of non- 
Hermiticity a — 7rj/(0)a, see Fig. |[ This fact signals on the gradual decorrelation of the real 
parts Re Zi of complex eigenvalues. It can be easily understood because of increasing possibility 
for eigenvalues to avoid one another along the Y — Im Z direction, making their projections on 
the real axis X to be more independent. 

In order to study the difference from the Hermitian case in more detail let us consider again 
the large behaviour. Then it is evident, that the number variance is only slightly modified by 
non-Hermiticity as long as a <C ia; . We therefore consider the case a 3> 1 when we expect essential 
differences from the Hermitian case. For doing this it is convenient to rewrite Eq.(^7|) as a sum of 
three contributions: 



.(3) 



{00 
1 2 fdk .3, 
]^^^^ (TTfcjexp 


00 

2 /"rffc . 2, ,x 
j —r- sm (ttk) exp 

TT^ J k 



00 

2 [dz 



ak 



(1 — z) sin^ (ttzLx) exp[— (az)^ 



(58) 



First of all we notice, that for large a the third contribution Sj^-* is always of the order 0(exp —a^) 
and can be neglected. 

The relative order of the first and second terms depends on the ratio L^/a. In a large domain 
1 <^ Lx ^ a, the second term S2^^ is much smaller than S2^^ This implies that the number 
variance grows like T,{Lx) — Lxf{Lx/a). We find it more transparent to rewrite the function f{u) 
in an equivalent form: 



Ittu V 



dte 



which can be obtained from Eq.(p8D after a simple transformation. 

For w = L.^/a <C 1 we have simply / ~ 1 and hence a linear growth of the number variance. For 
M » 1 we have / ~ (Tr^/^-u)^^. Thus, £2(^2;) slows down: J^2iLx) ~ 

(2) 

Only for exponentially large Lx such that In (Lx/a) ^ a the term produces a contribution 
comparable with S2^''. To make this fact evident we rewrite ^2^^ as: 



^2 — 



2 fdkln(^]^lsin'(^]e-''' 



(59) 
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For L/a ^ 1 we can neglect the oscillatory part of the integrand effectively substituting 1/2 
for Eq.(|59[). The resulting integral can be evaluated explicitly. Remembering that 



.(1) 



Ux>>a~ a/(27r^/^) we finally find: 

S2(ix > a) = 



t3/2 



^ In — - 4 



7 



where 7 is Euler's constant. This logarithmic growth of the number variance is reminiscent of that 
typical for real eigenvalues of the Hermitian matrices. 

Another important spectral characteristics which can be simply expressed in terms of the cluster 
function is the small-distance behavior of the nearest neighbour distance distribution |^,^,^ . We 
present the derivation of the corresponding relationship in Appendix. 

Substituting the expression Eqs.( p5|p4| ) for the mean density and the cluster function into 
Eq. ( [A^ ) one arrives after a simple algebra to the probability density to have one eigenvalue at the 
point Zq = X + iyoA and its closest neighbour at the distance \zi — zo\ — sA, A — {i'{X)N)^^ , 
such that s ^ 1: 

2' 



Zn 



9a{yo)^ga.{yo) - ( — 5a(2/o) 



d_ 

dyo 



exp ( -4| 



d0 exp 



(s^ cos^ ( 



2yos cos 6) 



(60) 



where ga is given by Eq. (p^) . 

First of all it is easy to see that in the limit a ^ 1 one has: Pa^i{Zo, s ^ 1) = ^{s/a^)^ in 
agreement with the cubic repulsion generic for strongly non- Hermitian random matrices ||3^ , [l^ , ^o|] . 
On the other hand one can satisfy oneself that in the limit a ^ we are back to the familiar GUE 
quadratic level repulsion: pa^o{ZQ,s ^ 1) cx S{yo)s'^. In general, the expression Eq.(|60|) describes 
a smooth crossover between the two regimes, although for any a 7^ the repulsion is always cubic 
for s 0. 

To this end, an interesting situation may occur when deviations from the Hermiticity are very 
weak: a <^ \f2 and 'observation points' Zq are situated sufficiently far from the real axis: 2|yo|/a ^ 
2-1/2. 

Under this condition the following three regions for the parameter s should be distinguished: 
i)- < TT^ ii)7r-r < - < 2^ and finally in) 2-1/2 ^ 2^ < ^ < a^i. 

In the regimes i) and ii) the term linear in cos Q in the exponent of Eq. (|60| ) dominates yielding 
the result of integration to be the modified Bessel function tt/q (^fl^)- In the regime iii) the term 
quadratic in cos6' dominates producing 27re-(*/''^^/o [(s/a)^] w (27ra/s)^^^. As the result, the 
distribution p(Zo, s) displays the following behaviour: 
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exp 



n2 



(61) 



with 50(2;) = 5a(2/)U=o- 

Unfortunately, the unusual power law p[s) oc s^/^ might be a very difficult one to detect numer- 
ically because of the low density of complex eigenvalues in the observation points reflected by the 
presence of the Gaussian factor in the expression Eq.(|6l|). 



IV. CONCLUSION 



In the present paper we addressed the issue of eigenvalue statistics of large weakly non-Hermitian 
matrices. The regime of weak non-Hermiticity is defined as that for which the imaginary part ImZ 
of a typical complex eigenvalue is of the same order as the mean eigenvalue separation A for the 
Hcrmitian counterpart. 

Exploiting a mapping to the non-linear cr— model we were able to show that there are three 
different "pure" classes of weakly non-Hermitian matrices: i) almost Hermitian with complex 
entries ii) almost symmetric with real entries and iii) complex symmetric ones. Within each of 
these classes the eigenvalue statistics is universal in a sense that it is the same irrespective of the 
particular distribution of matrix entries up to an appropriate rescaling. There are also crossover 
regimes between all three classes. 

Our demonstration of universality was done explicitly for the density of complex eigenvalues 
of matrices with independent entries. Within the non-linear it— model formalism one can easily 
provide a heuristic proof of such a universality for higher correlation functions as well as for 
" rotationally invariant" matrix ensembles, see [5^ ]. The above feature is a great advantage of the 
supersymmetry technique. 

A weak point of that method is a very complicated representation of the ensuing quantities. It 
seems, that the explicit evaluation of the higher correlation functions is beyond our reach at the 
moment, and even a calculation of the mean density requires a lot of effort, see [ p5|j27| ]. As a result, 
at present time the mean density is known explicitly only for the cases i) and ii). 

Fortunately, because of the mentioned universality another strategy can be pursued. Namely, 
one can concentrate on the particular case of matrices with independent, Gaussian distributed 
entries for which alternative analytical techniques might be available. Such a strategy turned out 
to be a success for the simplest case of complex almost-Hermitian matrices, where we found the 
problem to be an exactly solvable one by the method of orthogonal polynomials. This fact allowed 
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us to extract all the correlation functions in a mathematically rigorous way 

One might hope that combining the supersymmetric method and the method of orthogonal 
polynomials one will be able to elevate our understanding of properties of almost-Hermitian random 
matrices to the level typical for their Hermitian counterparts. 

/^From this point of view a detailed numerical investigation of different types of almost-Hermitian 
random matrices is highly desirable. Recently, an interesting work in this direction appeared 
motivated by the theory of chaotic scattering [Q . Weakly non- Hermitian matrices emerging in 
that theory are different from the matrices considered in the present paper because of the specific 
form of the skew Hermitian perturbation, see e.g. This fact makes impossible a quantitative 
comparison of our results with those obtained in [ [74[ . The qualitative fact of increase in number 
variance with increase in non-Hermiticity agrees well with our findings. Let us finally mention that 
the knowledge of the time-delay correlations [ pO| ]] allows one to make a plausible conjecture about 
the form of the number variance for the scattering systems with broken time-reversal symmetry. 
These results will be published elsewhere [ ||75| ]. 
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APPENDIX: NUMBER VARIANCE AND NEAREST NEIGHBOUR DISTANCE 

DISTRIBUTION. 

The number of eigenvalues in any domain A is expressed in terms of the eigenvalue density as 
in Eq.(|4l|). Then the variance T^2{A) = {n{A)^) — of the number n{A) is given by: 



S]2(A) = J d'Z^ J d'Z, [{p{Z,)p{Z2)) - {piZ,)){p{Z,))] 

A A 

= Jd^Z {p{Z)) - Jd^Zi jd^Z2 3^2(^1, Z2), 

A A A 

where we used the definition of the cluster function 3^2(^1, ^2)- 

We are interested in finding this variance for the domain A being a rectangular in the complex 
plane Z = X + lY: {) < X < Lx\ < Y < Ly. Moreover, we are going to consider the 

extension being comparable with the mean eigenvalue separation along the real axis: A = 
l/{h'{X)N). We know that on such a scale the mean eigenvalue density is independent of X and 
can be replaced by its value i^(0) at X = 0, whereas the cluster function depends on fl ~ Xi — X2 
rather than on Xi and X2 separately. Using these facts, we obtain: 

Ly Lx Ly Ly 

^2{L^;Ly) = L, J dY {p{0,Y))-2jduj {L^-n) J dYi J dY2 y2m,Yi,Y2) (Al) 
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It is convenient to introduce the spectral form-factor B{K, Qi, Q2) by the Fourier transform: 

C30 00 00 

y2i\niYuY2) ^ J dK JdQi y^ga i?(if,Qi,Q2)e-2"("™'5^+^=^«=) (A2) 



-00 —00 



The number variance can be expressed in terms of the spectral form-factor, Eq.(A2), as: 

L/Tj 



Ly) =L, J dY {p{0, ¥)) - (A3) 

— Ly 

oc 00 00 

2 f f sin27rQisin27rQ2 f .r^ ^''^^^ ^^^^ ut n it r, i t \ 
— y dQi j dQ2 J dK — B{K,Qi/Ly,Q2/Ly). 

— 00 ~CKj 

In particular, for the strip < X < L; —00 < y < 00 the number variance is given by a rather 
simple expression: 

°° ■ 2 

00) = LNiy{0) [dK 0, 0) (A4) 



In the main text of the paper we use the variables 1/1^2 — NYi_2 and w — Nfl. Correspondingly, 
the form- factor B{K, Qi, Q2) is related to h{k, gi, 52), Eq.(^) as 

B{K, 0, 0) = -^KK/N, 0, 0) = iVe— '"'(^/^)' (zy(0) - \K/N\) (j/(0) - \K/N\) 



Substituting this expression into Eq.(A4) and measuring the length L in the units of A = 
l/(A^i'(0)) as L = ALx we find after simple manipulations the eq.(^7[). 

Let us now derive the relation between the cluster function and the nearest neighbour distance 
distribution ^(^0,5"), see also pO| ]. 

We define the quantity p{Zo,S) as the probability density of the following event: i) There is 
exactly one eigenvalue at the point Z = Zq of the complex plane, ii) Simultaneously, there is 
exactly one eigenvalue on the circumference of the circle \Z — Zq\ = S iii) All other eigenvalues Zi 
are out of that circle: — Zo| > S. 

00 

As a consequence, the normalization condition is: j d^Z^ J dS p{Zq, S) — 1. In particular, for 



Hermitian matrices with real eigenvalues one has the relation: p{Zq,S) = S(lmZQ)v{XQ)px{S), 
with px{S) being the conventional "nearest neighbour spacing" distribution at the point X of the 
real axis ||]. 

Using the above definition one easily finds the relation: p{Zo, S) = —jgH{Zo, S), where H{Zo, S) 
has the meaning of the probability density to have one eigenvalue at Z = Zq and no other eigenval- 
ues inside the disk D:\Z — Zo\ < S. The latter quantity is related to the joint probability density 
of complex eigenvalues as 

N 

H{Zo,S) = ^Yl J d'^Zi... j d^ ZnVn{Zi, ...,Zn)5^'-\z, - Zo)\{[l - xd{Z,)] (A5) 
1=1 
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where xoiZj) is the characteristic function of the disk equal to unity for points Zj inside the disc 
and zero otherwise. 

We are interested in finding the leading small-S' behaviour for the function p{Zo, S). For this 
one expands 

n [1 - Xd{Z,)] = l-J2 Xd{Z,) + E E XD{Z,)xD{Zk) - ... 

and notices that each factor xoiZj) produces upon integration extra factor proportional to the 
area of the disc. Therefore, to the lowest nontrivial order in S one can restrict oneself by the first 
two terms in the expansion and write: 

1 



H{Zo,S) ^ 



{piZ)p{Zo)) - / d'ZxD{Z)R2{Zo,Z) 



where we used the definitions of the mean eigenvalue density and the spectral correlation function, 

S 2iT 

see Eq.®. At last, exploiting that / (PZxdF{Z) = J drr J d0F {Zq + re''^) one finally finds 



after differentiation over S: 

p{Zo,S)^^ J d9[{p(Zo)){p{Zo + Se''))~y2{Za,Zo + Se'')] (A6) 



where we used the definition of the cluster function, Eq. (^3|) . 

In the regime of weak non-Hermiticity this formula is valid as long as the parameter S is small 
in comparison with a typical separation between real eigenvalues of the Hermitian counterpart: 
S'< A~ 1/iV. 
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FIG. 1. The cluster function Y2, defined by Eq. (g5|), against cj for the parameter values yi = y2 = 
and a — (left solid line, the GUE case), a = 0.5 (dashed line), a = 1.0 (dot-dashed line), a — 1.5 (dotted 
line), and a — 2.0 (the right solid line, well approximates the Ginibre case) 
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FIG. 2. The number variance E2, defined by Eq. against Lx for the parameter values a = (sohd 
fine, the GUE case), a — 1.5 (dashed hne), and a = 2.0 (dot-dashed line) 



